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We have analyzed a sample of nearby cool and warm infrared (IR) galax- 
ies using photometric and structural parameters. The set of measures include 
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surface brightness as well as radio, near-infrared, and optical sizes. In a given 



are found approximately la below and above the mean color in the far-infrared 
C — Ltir diagram. We find that galaxy radio surface brightness is well cor- 
related with color whereas size is less well correlated with color. Our analysis 
indicates that IR galaxies that are dominated by cool dust are large, massive spi- 
rals that are not strongly interacting or merging and presumably the ones with 
the least active star formation. Dust in these cool objects is less centrally con- 
centrated than in the more typical luminous and ultra-luminous IR galaxies that 
are dominated by warm dust. Our study also shows that low luminosity early 
type unbarred and transitional spirals are responsible for the large scatter in the 
C — Ltir diagram. Among highly luminous galaxies, late type unbarred spirals 
are predominately warm, and early type unbarred and barred are systematically 
cooler. We highlight the significance of C — Ltir diagram in terms of local and 
high redshifts sub-millimeter galaxies. 

Subject headings: galaxies: general - galaxies: spirals - galaxies: active - galaxies: 
starburst - galaxies: structure 
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1. Introduction 

Observations by the Infrared Astronomical Satellite (IRAS) led to the identification of 
very luminous and ultra-luminous infrared galaxies, commonly known as the LIRGs and 
ULIRGs, respectively. These galaxies have enormous far-infrared (FIR) luminosities, 10 11 < 
L FIR /L Q < 10 12 for LIRGs and L FIR > 10 12 L for ULIRGs, and emit the bulk of their 
energy at the infrared wavelengths (Soifer et al. 1984, 1986, 1987; Sanders et al. 1988). 
Subsequent multi- wavelengths studies reveal that in the local universe (z < 0.1) ULIRGs 
have comparable total luminosities (Ltir) and higher space densities than those of optical 
quasars (Sanders et al. 1988; Kim & Sanders 1998). 

Understanding the nature and origin of energy sources in the LIRGs and ULIRGs has 
been the subject of much debate. Studies indicate that in the vast majority of these objects 
power comes from dust which is heated to various temperatures by various thermal and non- 
thermal processes such as ongoing (steady-state) star formation, intense starburst phase, 
and/or synchrotron radiation from the supernovae explosions (Lawrence et al. 1986; Genzel 
et al. 1998; Lutz et al. 1998). These systems may also be powered by active galactic nuclei 
(AGN; Lonsdale, Smith, & Lonsdale 1993). The energy distributions suggest that most of 
the observed FIR emission from luminous disk galaxies is due to the thermal radiation from 
warm dust heated by hot stars embedded in HII regions and molecular clouds heated directly 
by young OB stars; and cool dust from the general interstellar radiation field (Helou 1986; 
Presson & Helou 1987; Barvainis, Antonucci, & Coleman 1992). Dust could also be heated 
by shocks in the interstellar medium during the collision or interaction of galaxies (Harwit 
et al. 1987). 

Dust heating mechanisms can be traced indirectly as each of these should correspond to 
particular spatial distributions of the FIR light. For example, the FIR emission originated 
by galaxies dominated by an AGN would appear compact and unresolved. On the other 
hand, starburst heated dust should have about the same scale size as the burst itself because 
young stars are well mixed with the gas from which stars are forming. Dust heated by 
the non-ionizing photons from cold, population I stars, might be expected to follow the 
smoothed distribution of older giants in the galaxy (Zink et al. 2000). A determination of 
the total energy output such as luminosity, the amount of dust present in the interstellar 
medium and its spatial distribution within the galaxies as well as its relationship to other 
basic components, such as the atomic and molecular gas, the stars, etc, is essential to studies 
connecting galactic structure and nature of the IR energy sources (Carico et al. 1990; Sopp 
& Alexander 1991, 1992; Andreani & Franceschini 1992, 1996; Alton et al. 1998; Domingue 
et al. 1999; Siebenmorgen, Kriigel & Chini 1999; Haas et al. 2000; Stickel et al. 2000; 
Trewhella et al. 2000; Zink et al. 2000). 
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The FIR color C, defined by the IRAS 60 -to- 100 fim flux density ratio, an indicator 
of characteristic dust temperature, is a diagnostic of the typical heating conditions in the 
interstellar medium of a galaxy (Bothun, Lonsdale, & Rice 1989; Soifer & Neugebauer 1991). 
It has been shown that low-redshift IRA S galaxies exhibit positive correlation of C with Lpm 
such that more luminous galaxies tend to be warmer as compared to less luminous galaxies 
(Soifer et al. 1987; Dale et al. 2001; see Fig.l in this study). The color changes systematically 
over roughly three orders of magnitude in luminosity with large scatter. As a result there 
is a substantial number of highly luminous but cool galaxies as well as low luminosity yet 
warm/hot galaxies (Dale et al. 2001; Chapman et al. 2003). In this study we attempt to 
understand cool and warm IR luminous sources in terms of the geometric distribution of 
dust and galaxy morphology. In a given FIR luminosity range, these galaxies are considered 
in the context where cool galaxies are \a below and warm galaxies are la above the mean 
color at a given luminosity. We are interested in differences between cool and warm galaxies 
in the same luminosity range that might suggest differences in the nature of the heating 
sources. 

In an interesting study Young (1999) showed that star formation efficiency (SFE), de- 
fined by the ratio of Ltir to molecular hydrogen mass Mh 2 , anti-correlated with disk (optical) 
sizes. The trend is prevalent in galaxies of various Hubble types and environments. It has 
been attributed to the shear present in disk galaxies. The molecular clouds in larger disks 
would experience increased turbulence which would reduce the efficiency of star formation 
(Young 1999). Galaxies in Young's sample span a broad range in FIR color comprising 
normal star forming galaxies (SFGs, L FIR / L Q < 10 11 ) as well as LIRGs and ULIRGs. In 
Young's study, however, it remains unexplored how cool and warm systems behave in the 
visual as well as in the longer wavelengths. Since the IR luminous galaxies contain very 
large amount of dust, the interpretation of any observation at optical or shorter infrared 
wavelengths gets complicated by extinction. Therefore, to probe the origin of large IR lumi- 
nosity it is necessary to have the knowledge of the spatial extents of the emitting regions in 
observations unaffected by dust extinction. 

The paper is organized as follows: we describe the data in §2 and parameter error 
estimation in §3. We present our results in §4 and discuss the implications of our results in 
§5. The conclusions are given in §6. 

2. Data 

IRAS resolved only the largest and nearest galaxies because of its comparatively large 
beam size of ~ 2 -4 at 60 and 100 jim. Therefore, the FIR brightness distributions of most 
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of the sources detected by the IRAS are unknown. In contrast to the rather large scatter in 
optical-FIR correlations, a tight correlation of FIR with radio continuum total flux densities 
for infrared-selected galaxies appear to hold locally within individual galaxies (Helou, Soifer 
& Rowan-Robinson 1985; Beck & Golla 1988; Bicay, Helou, & Condon 1988; Murphy et al. 
2006). In nearby galaxies where (IRAS) FIR sizes have been measured directly, radio sizes 
are known to match, or to be somewhat larger than the size scale in FIR (Bicay & Helou 
1990). Besides, FIR and radio continuum brightness distributions of these galaxies show 
remarkable similarity (Marsh & Helou 1995). The agreement in FIR and radio sizes, albeit 
derived from a smaller sample, suggest that high resolution radio maps may be taken as a 
good substitute, or at least as an upper limit for the size of the unobtainable FIR maps. We, 
therefore, use 1.49 GHz radio size as a proxy for the size of FIR emission. 

Wang & Helou (1992, hereafter WH) studied the compactness of FIR bright galaxies 
and constructed a list of 330 galaxies from the IRAS Bright Galaxy Sample (BGS; Soifer 
et al. 1986, 1987, 1989). WH constructed a flux limited sample of extragalactic objects 
brighter than 5.24 Jy at 60 fim, covering the entire sky surveyed by the IRAS at galactic 
latitude |6| > 30°. For radio data (rest frame flux density, size) WH used the atlas of the 1.49 
GHz radio maps based on VLA observations of IR bright galaxies that had been compiled by 
Condon, Anderson & Helou (1990, hereafter CAH). After removing AGN dominated galaxies 
identified as highly compact radio sources, the final list of WH contained 218 galaxies. Our 
sample is based on this list of 218 galaxies. We use the same radio data, however, we take the 
rest frame IRAS flux densities, luminosity, and estimated distance of member galaxies from 
the Revised Bright Galaxy Sample (RBGS) since it provides the best available reference for 
accurate IRAS fluxes and IR luminosities of galaxies in the local universe (see Sanders et 
al. 2003 for details on RGBS). 

Particular attention is given to exclude AGN from the sample. Although the list of WH 
was carefully compiled to get rid of the radio monsters, we have checked the entire sample 
following the prescription of de Grijp et al. (1985) using mid and far infrared color - color 
correlations. We have found 5 galaxies that fall into the region of the color - color plot which 
is mostly occupied by AGN like sources (de Grijp et al. 1985). We have removed them from 
the list. The radio fluxes of the remaining 213 galaxies correlate well with the IRAS fluxes, 
and all galaxies fall within ~ 1.5er of the mean q- value in the L FIR — q correlation indicating 
absence of radio excess objects (q parameter is defined in Condon, Anderson & Helou 1991). 
We have also made visual inspection of the observed spectral energy distributions of member 
galaxies. From this analysis we conclude that our final sample of 213 sources are star forming 
disk galaxies, although there may be AGNs present that are not energetically dominant in 
the far- infrared. 
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The number of SFGs, LIRGs, and ULIRGs in our sample are 154, 52, and 7, respectively. 
These three classes of galaxies have the following redshift distributions: 0.0 < z < 0.016, 
0.012 < z < 0.051, and 0.018 < z < 0.082 with the median redshifts 0.006,0.023, and 0.055. 
The effective color temperature of the integrated dust emission in the sample galaxies ranges 
in between ~ 25K — 40K, assuming emissivity index, (3 = 2 (Dunne et al. 2000; Dunne & 
Eales 2001). The ULIRGs show signs of interaction, e.g., presence of neighbor and accretion, 
e.g., tidal tails and disturbed outer envelopes. Majority of LIRGs appear as single, isolated 
systems while 30% of them appear disturbed. At the lower end of the luminosity range, the 
LIRGs are luminous single isolated galaxies, with features resembling to the low luminosity 
SFGs. A minority of SFGs (15%) has disturbed outer envelope or nearby companion. 

The NASA/IPAC Extragalactic Database (NED) and the Lyon-Meudon Extragalactic 
Database (LEDA) has been used to obtain Galactic extinction corrected S-band luminosity 
L B and 25 mag arcsec -2 linear diameter D B - The near-infrared (NIR) K s band (2.17/im) 
magnitude of our sample is 5.0 < K s < 10.6. For bright galaxies Jarrett et al. (2000) have 
recommended K s band 20 mag arcsec" 2 diameter Dnir as the most reliable and robust for 
galaxy photometry. We obtain D^ir from Two Micron All Sky Survey catalog (Jarrett 
2000) using a 5" search radius. We did not find any reliable 5-band size for CGCG247- 
020, IRASF 08339+6517, and IRASF 12132+5313, and NIR size for NGC 5256, NGC 5331, 
and MCG+07-23-019. All these galaxies belong to the LIRGs sub-sample. No correction 
was applied to the luminosities and diameters for internal extinction or inclination of the 
host galaxies. Following RBGS we adopt Q M = 0.3, JIa = 0.7, and H = 75 km sec -1 
Mpc -1 . These values are slightly different than recent estimates from Wilkinson Microwave 
Anisotropy Probe (WMAP; Spergel et al. 2003). Table 1 includes a partial list of low 
luminosity star forming galaxies showing only cool and warm sources. In the table galaxies 
are sorted in ascending order of FIR color, i.e. galaxy numbered 1 has the lowest FIR color, 
for all luminosity classes [See the electronic edition of the Journal for the complete list of 
galaxies] . 



3. Uncertainty in Parameters 

We use Ltir as a photometric measure, instead of Lpm, since the former is based on 
the all four fluxes measured by the IRAS and thus taking contributions from almost the 
entire IR range (Sanders & Mirabel 1996). The median uncertainties associated with the 
IRAS flux densities are ~ 5% and ~ 3% at 12 and 25 microns, respectively, and ~ 1% in 
both far-infrared bands. This leads to a 2% uncertainty in the FIR color of a galaxy which 
is much smaller than the spread in mean color in a given luminosity. The total IR flux of 
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a galaxy is a weighted sum of four IRAS flux densities. With the corresponding weighting 
factors given in Sanders & Mirabel (1996) the median uncertainty in the total IR flux of a 
galaxy is ~ 1%. 

Distances for the galaxies are taken from the RBGS. Most of the RBGS distance esti- 
mates come from redshift measurements, application of the Hubble law and correction for 
the Mould et al. (2000) flow model. However, some are primary (P) or secondary (S) dis- 
tance estimates as flagged in the RBGS Table 1 which do not come from the Hubble flow 
and cosmic attractor model. Our sample contains 48 galaxies of such distance estimates, 
all of which belong to the SFGs sub-sample. A significant uncertainty to a galaxy size and 
Ltir would result from an uncertainty in the distance estimate. Surface brightness being 
the distance independent measure would not be affected. The observed heliocentric radial 
velocities of sample galaxies are measured ~ 1 — 2% accuracy (obtained from NED). The 
uncertainties associated with the model parameters such as motions in the Local group, 
Virgo centric infall, Great Attractor infall etc. are high (> 10%; Mould et al. 2000). The 
uncertainty associated with other distance estimates is of similar magnitude. A galaxy op- 
tical angular size has less than 10% uncertainty, depending on the flattening of the disk, as 
given in the Third Reference Catalog (de Vaucouleures et al. 1991; RC3). We assume a 
similar uncertainty in the NIR and radio bands (we note that it varies in different bands). 
Combining all of the preceding, we assign ~ 15% and ~ 20% uncertainty, respectively, to 
the estimate of physical size and Ltir of a galaxy. 

The uncertainty in 1.49 GHz surface brightness is due to the measurement errors in 
1.49 GHz flux density and angular size. According to CAH (1990), the contributions from 
rms confusion error and calibration error to the 1.49 GHz flux density are smaller than 
the noise error. The rms radio map noise a n is between 0.1 and 0.2 mJy per beam solid 
angle Vl b ~ 9\ eam (Condon 1987; CAH 1990). The noise contribution to the flux density is 
a s ~ a^Vl/VLb) 1 / 2 in mJy, where f2 ~ 9^ x 9^ is the solid angle covered by the source and 
®M-> ®m are the deconvolved major (M) and minor (m) axes of a radio image. Radio maps in 
CAH (1990) ranges in FWHM angular resolution from 1.5 ' to 60 " depending on the source's 
apparent size. Using a n = 0.2 mJy for each galaxy gives an uncertainty of ~ 2% in 1.49 GHz 
flux density. Uncertainty in radio size varies with beam resolutions but it is within 5% in all 
cases. Taking this upper limit we find ~ 10% uncertainty in the surface brightness coming 
from 1.49 GHz flux density and angular size. 
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4. Results 

Figure 1 shows the FIR color - IR luminosity diagram for IRAS Seo^m > 1.2 Jy sample 
of ~ 4700 sources (Fisher et al. 1995) in panel a and our sample of 213 galaxies with 
Sm^m > 5.4 Jy in panel b. The former sample is shown after removing 330 spurious cold 
luminous galaxies as pointed out by Chapman et al. (2003). In both panels horizontal bar 
and small cross represent median, and mean color with la error bar. These panels highlight 
the fact that galaxy FIR color and total IR luminosity follow a broad correlation with large 
spread in different flux limited samples. In this figure small triangles and dots represent SFGs 
and LIRGs, respectively (Fig. lb). For these two classes of galaxies solid symbols represent 
cool sources (around or below la from the mean C) and open symbols represent the warm 
sources (around or above la from the mean C. Open stars are used for the ULIRGs. Because 
of their limited number we are unable to separate these ultra-luminous galaxies into cool 
and warm categories. In this panel large crosses represent local sub-milimeter galaxies (see 
below for a discussion of these galaxies). 

Figure 2 shows color as a function of SBiA9GHz as well as galaxy size in radio, NIR, 
and optical bands. The choice of these parameters are based on conventional practice which 
suggests that radio size is a stand-in for FIR emission size, NIR size represents emission 
from old stellar population, and optical size accounts for the emission from a mixture of 
both young and old stars but complicated by dust. We have used the major axis to estimate 
the area of each object to reduce the effects of highly uncertain inclinations. The radio 
surface brightness is calculated as Sia9GHz/{0m/2) 2 - 

Cool and warm galaxies separate clearly in the surface brightness vs. color diagram 
(Fig. 2a). For each luminosity class, warm colors correlate with high surface brightness. 
Panel b indicates that galaxies that are dominated by cool dust are large, massive spirals 
that are not strongly interacting or merging and presumably the ones with the least active 
star formation. Dust in these cool objects is also, on average, less centrally concentrated 
than in the more typical LIRGs and ULIRGs that are dominated by warm dust. There is a 
clear sequence in the progression of color with radio size: color increases systematically from 
normal/isolated disks toward merger/disturbed disks. The latter type of disks are shown by 
asterisk symbols. Panels a and b shows that color is strongly correlated with both surface 
brightness and galaxy size at 1.49 GHz (see Table 2 for correlation statistics). Note that 
the uncertainty in color or brightness is statistical in nature and do not include calibration 
error. However, it will not affect our results since incorporating calibration error will simply 
shift these physical parameters systematically. 

The statistical trend in Fig. 2b and Fig.2c suggests that luminous cool galaxies, where 
older stars have greater spatial distribution, show a tendency to have cool dust distributed 
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over larger volume. This may be interpreted, as an indirect but interesting support, that FIR 
light from cool galaxies are reprocessed emission of photons coming predominantly from the 
old population of stars. This is in accord with the expectation that dust heated by the non- 
ionizing photons from a cooler population of stars might follow the smoothed distribution 
of older stars in the galaxy. In spite of a large scatter, we can see a moderate trend where 
warm sources have relatively smaller regions of NIR emission than their cool counterparts. 
Most of the ULIRGs (4/7) have larger NIR and 5-band diameters than found among warm 
(C > -0.2) LIRGs (Fig.2c and 2d). It is interesting that ULIRGs, which all have very 
warm dust (C > —0.1), have NIR and B diameters comparable to the subset of LIRGs with 
relatively cool dust temperatures (C < —0.22). The data indicate that the ULIRGs have 
relatively large disks (e.g., they involve mergers of massive galaxies) with compact cores 
that dominate the FIR emission. It is also clear from Fig. 2c that warm LIRGs (large open 
circles) have systematically smaller NIR diameters than cool LIRGs (large solid circles). We 
conclude that some galaxies with relatively cool dust temperatures are LIRGs (rather than 
lower luminosity SFGs) because, despite having a lower average SFR per unit area, they 
have on average larger total surface areas than warm LIRGs. 

The Optical size, representing the spatial distribution of stellar emission, shows no 
trend with the color when the whole sample is considered (Fig.2d). However it improves 
for SFGs+LIRGs, after removing ULIRGs from the sample (see Table 2 for correlation 
statistics). The scatter is large as compared to Fig.2b and Fig. 2c. All ULIRGs and a 
substantial number of warm SFGs and LIRGs are almost similar in optical size compared 
to their cool counterparts. All of these warm sources show signs of interaction such as tidal 
tail, disturbed outer region or presence of neighbors. 

With the exception of color and SBi A9GHz all of the physical parameters depend on the 
distance. Figure 2 demonstrates that color is tightly correlated with distance independent 
measures. On the other hand correlation between color and galaxy size is weak. There are 
at least three factors that can contribute to the broader distribution in color-size correla- 
tions (see Table 2 for correlation statistics). First, there may be large error in the distance 
measurement in spite of correction for non-Hubble flows since correlation between color vs. 
NIR or B-band surface brightness is relatively stronger compared to the respective sizes. 
Second, optical diameter D B of interacting/merging galaxies will overestimate the actual 
size because of the inherent difficulty associated with size estimates of these sources. Re- 
moving the ULIRGs from the sample makes the color-size trend stronger. Third, varying 
scatters in color vs. size plots reflect methodological differences in estimating angular size in 
different wavelengths. The scatter would have been reduced if galaxy sizes were estimated 
in all wavelengths in a systematic manner, e.g., at the same brightness level. 
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Figure 3 shows color as a function of galaxy (optical) morphology. We use RC3 classifi- 
cation and divide galaxies into three broad categories: unbarred (A), transitional (AB), and 
barred (B). In each morphology bin, we take S0?-S0a-Sa-Sb-Sbc sources as early type, and 
S?-Sc-Scd-Sd-I-Pec sources as late type. Early and late type division is made to be consistent 
with earlier studies that have found fundamental differences in the properties of early and 
late type galaxies, especially, barred galaxies (Combes & Elmagreen 1993; Ho, Filippenko & 
Sargent 1997; Sakamato et al. 1999; Sheth et al. 2005). 

We did not find any classification for 11 LIRGs and it is also uncertain for the seven 
ULIRGs [See Table 1 in the electronic edition of the Journal for these galaxies]. As a result, 
Fig. 3 shows 154 SFGs and 41 LIRGs divided into SA, SAB, and SB type sources. This 
figure illustrates that high luminosity LIRGs are mostly unbarred spirals (Fig.3a, b). The 
transitional and barred spirals, in general, fall into the class of low luminosity SFGs (Fig. 
3c,d and e,f). The essence of this figure are the followings: for SFGs, late type spirals stay 
close to the mean color (Fig. 3d and f). Early type spirals, on the other hand, are responsible 
for the scatter in color-luminosity diagram (Fig. 3a, c, and e). Interestingly, cool sources 
consist of early type unbarred and transitional galaxies whereas warm sources consist of 
early types of all three categories. In other words, low luminosity barred galaxies, in general, 
show a tendency to have higher color temperature. For LIRGs, late type unbarred galaxies 
are predominately warm (Fig. 3b). This is not unexpected since late type galaxies are highly 
irregular systems with higher star formation rate which results in warmer color. Early type 
LIRGs are systematically cooler (Fig. 3a, c, and e). Both transitional and barred spirals are 
extremely rare in this luminosity class. 

In a recent study, Sheth et al. (2005) showed that late types barred spirals are less 
centrally concentrated than early types and a significant subset of early type barred spirals 
have little or no gas within the bar region. This observation has been explained as a result of 
higher mass accretion rates in the past in early type barred spirals where the large amount 
of gas driven inward by the bar have already been converted into stars. They suggested that 
these galaxies are in the post-starburst phase. The IR luminosity of Sheth et al. sample of 44 
galaxies span the range, 10 9 < Ltir/L & < 10 11 , which we call star forming galaxies (SFGs) 
in our sample. These galaxies fall in the class of SFGs in our sample. The observation that 
early type barred spirals are in the post-starburst phase gives a plausible explanation for 
these galaxies to be low luminosity star forming systems. 

Note that the RC3 classifications used here are based on optical images that are strongly 
affected by dust obscuration in our FIR-selected galaxy sample. More recent studies at NIR 
wavelengths indicate the presence of stellar bars in a large fraction of disk galaxies that appear 
unbarred at optical wavelengths (Eskridge et al. 2000). In addition, when one extends the 
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bar detection sensitivity to a low relative amplitude of ~ 3% in the K band, nearly 90% of 
a sample of optically unbarred (SA) spirals contain stellar bars (Grosbol, Patsis, & Pompei 
2004). Therefore, the fractions of transition (SAB) and barred (SB) galaxies is likely to 
be much higher than is indicated by the RC3 classifications, and the optical bar fractions 
may be telling us more about the relative extinction by dust at optical wavelengths than 
the intrinsic stellar distributions. An improved understanding of this issue requires detailed 
analysis of near-IR images for a large sample of LIRGs and ULIRGs. 



5. Discussions 

It is necessary to address whether the large scatter in the FIR color-luminosity diagrams 
is real (intrinsic to the physics of the objects), or an artifact of the observations. Figure 1 
shows the data from the IRAS sample of Fisher et al. (S^m > 1-2 Jy; panel a) and data 
from the RBGS (S 60fim > 1.2 Jy; panel b). The main objective of the RBGS was to produce 
more accurate fluxes for the many large, nearby galaxies resolved by IRAS, by recovering 
extended emission that was not represented in the (underestimated) flux densities reported 
for such objects in the IRAS Point Source Catalog (PSC) and Faint Source Catalog (FSC). 
The RBGS also added some objects missed in the previous BGS compilations. This problem 
is minimal for the fainter objects in the 1.2 Jy sample, which on average are more distant 
and therefore smaller than the RBGS objects, with little or no flux missed in the published 
PSC/FSC flux densities. The main conclusion we can draw from comparing Fig. la with 
Fig. lb is that the scatter in the 5.24 Jy RBGS is consistent with the scatter in the Fisher 
et al. 1.2 Jy sample. The presence of more true statistical outliers (> 3a in each luminosity 
bin) in the 1.2 Jy sample compared to the 5.24 Jy sample is expected from basic sampling 
statistics, as is the extension to higher L TIR values due to sampling of a much larger volume 
of space in the 1.2 Jy sample. As stated in previous sections, the uncertainty in the Ltir 
values are dominated by the distance determinations (~ 15%), and the uncertainty in C is 
relatively small (See the representative error bar in Fig. la). Therefore, we conclude that 
the dispersion in the C — Ltir digram is intrinsic and related to the diversity of physical 
conditions of the ISM in the galaxies. 

To understand the physical origin of the color-luminosity diagram it is necessary to 
have detailed knowledge of these cool and warm galaxies. Questions can be raised such as, 
in a given luminosity range, what are the possible (external or internal) mechanisms that 
would cause some galaxies to have cool color temperature compared to the majority of the 
galaxies ? what are the basic connections between the photometric and structural properties 
of galaxies in this diagram ? We proceed in light of above queries and find that dust in the 
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cool objects is less centrally concentrated than in the more typical LIRGs and ULIRGs that 
are dominated by more centrally concentrated warm dust. We also find that the optical disks 
of IR cool galaxies show a tendency to be more extended than those of the warmer ones. The 
color-size trend is relatively stronger at longer wavelengths. The extended emission in cool 
sources may well indicate heating by old stars, but it could also indicate simply wide-spread 
small star forming regions scattered in a large dusty disk. Distinguishing between these 
aspects is beyond the scope of this study. 

The color-luminosity trend shown in Figure 1 is observed in the local universe (z < 
0.082). It had been demonstrated that the trend does not change out to z ~ 1 (Chapman 
et al. 2003) meaning that the dispersion of color in each luminosity bin does not vary 
significantly with redshift. However, given the presence of higher luminosity objects such 
as hyper-luminous galaxies with Ltir/L q > 10 13 in high-redshift samples that cover much 
larger volumes of space and earlier look-back times than local samples, from the present 
data we can only make such inferences for objects with L TIR /L & < 10 12 5 . In this respect, 
this correlation, therefore, reveals something quite profound: in absence of any structural 
information of high z luminous galaxies we can get a rough idea about the extent of the 
disks if we simply know the FIR colors and luminosities of these galaxies. In spite of a 
large scatter, the trend suggests that cool galaxies are generally highly luminous because 
of extended disks. The warm galaxies may or may not have extended disks depending on 
galaxy environment and/or internal processes. 

Current sub-millimeter /millimeter (sub-mm/mm) surveys have discovered a new pop- 
ulation of z > 1 sub-mm galaxies, contributing up to 50% of the extragalactic background 
light (Smail, Ivison, & Blain 1997; Barger et al. 1998; Hughes et al. 1998; Borys et al. 
2003). There are two possible explanations of these galaxies: they are either proto-spheroids 
with high star formation rate > 10 3 M yr _1 (Dunlop 2001) or disks dominated in the 
sub-millimeter by infrared cirrus heated by their interstellar radiation field, rather than in- 
tense star formation (Rowan-Robinson 2001; Efsthathiou & Rowan-Robinson 2003, hereafter 
ERR; King & Rowan- Robinson 2003). ERR have modeled cirrus emission for a selection of 
local (z < 0.02) and high-z (> 1) sub-mm galaxies and found excellent agreement in all 
cases. Their work suggests that sub-mm count and background can be understood in terms 
of cirrus-like emission (with effective dust temperature < 30 K) rather than invoking dusty 
ultraluminous starburst galaxies, e.g. Arp 220 in our local universe, to be the representative 
of distant sub-mm galaxies. The significance of ERR work is that if the cool component 
dominates the observed sub-mm fluxes of these high-z sub-mm/mm wave survey sources the 
emission should be extended rather then compact (i.e. centrally condensed). We attempt to 
explore this aspect of sub-mm galaxies in the context of our study by analyzing the 1.49 GHz 
emission of the local sample of ERR since only these galaxies are available in the RBGS. 
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The local sample of ERR includes UGC 0903, NGC 0958, NGC 1667, NGC 2990, UGC 
5376, NGC 5962, NGC 6181. We remove NGC 2990 from the ERR sample because it has 
Sm^m < 5.24 Jy placing it below the flux limit of the RBGS. Except NGC 1667 and UGC 
5376, the other four galaxies are present in our sample. The IR and radio information of 
these two galaxies are obtained respectively from the RGBS and the NRAO VLA Sky Survey 
catalog (NVSS; Condon et al. 1998). The FIR color and total IR luminosities of these six 
galaxies (shown by large crosses) are similar to normal SFGs and cool LIRGs in our sample 
(Fig. lb). We find that these galaxies indeed have large radio (FIR) disks. The 1.49 GHz 
radio source diameters of these galaxies range ~ 4 to 20 kpc, which is comparable to cool 
LIRGs and even normal SFGs (those not defined in this study as unusually cool or warm). 
This is in stark contrast to the warm SFGs and LIRGs that have radio sizes less than ~ 3 
kpc (See Fig.2b). 

The FIR cool high-z sub-mm galaxies of ERR sample carry a few xlO 13 L , several 
hundred times more luminosity than these local examples. Assuming that these galaxies 
are scaled-up versions of local cool examples, where disk size increases at constant surface 
brightness to account for the increased luminosity, we find that their radio (and FIR) disks 
must range up to 100 kpc or more in size. While rare, local galaxies with similar sizes 
do exist, as in the case of low surface-brightness, Hi-rich galaxies (Sprayberry et al. 1995; 
Bothun, Impey, & McGaugh 1997; Matthews, van Driel, & Monnier-Ragaign 2001). At 
earlier epochs the cool sub-mm galaxies may be single large disks able to regulate their star 
formation at the modest level required for these luminosities and FIR colors. But what do 
we know from observations about the extent of FIR/Radio emission in these distant sub-mm 
sources, especially the FIR cool ones ? 

For the assumed cosmology in this study, the linear disk size of 100 kpc corresponds 
to an angular size of > 7" for z > 1 sub-mm galaxies. On the other hand, the reported 
sizes of sub-mm galaxies in the redshift range 1 < -2 < 4 are significantly smaller, < 
(Ivison et al. 2002; Iono et al. 2006; Tacconi et al. 2006). This suggests that the objects 
measured are best understood as scaled-up versions of local FIR warm ULIRGs rather than 
cool extended disks (Tacconi et al. 2006), since nothing is known about their FIR colors. 
While cold ULIRGs had been detected in the redshift range 0.4 < z < 1, very little is known 
about their morphology or geometry (Chapman et al. 2002). However, the recently reported 
optical properties of BzK-15504 (Genzel et al. 2006) make it a very likely FIR cool luminous 
galaxy. At z — 2.38, its disk extends out to ~ 10 kpc, and has elevated distributed star 
formation throughout in addition to an accreting central black hole. The estimated total 
star formation rate for this galaxy is 140^gQ° M yr _1 (Genzel et al. 2006) suggesting that 
it is an ultra-luminous infrared galaxy with Ltir > 10 12 L . 
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The relative frequency of occurrence of FIR cool objects as a function of redshift is still 
unknown, and the detection efficiency will decrease as the source extent increases. Tacconi 
et al. have measured 8 objects, all of which were compact. Their results only place a limit 
on the effective dust temperature (39 ±3 K for j3 = 1.5) of sub- mm galaxies rather than the 
frequency of cool objects. As precise redshifts and more accurate rest frame FIR fluxes of 
high-z galaxies become available, placing these sub-mm/mm galaxies in the context of the 
FIR color-luminosity diagram will help shed light on their structure. Combining that with 
kinematic information may constrain further the nature and evolving path of these galaxies. 

Warm LIRGs and ULIRGs are usually associated with intense bursts of star formation. 
These galaxies are known to have higher star formation rates (SFR) and star formation 
efficiency (SFE) than lower luminosity galaxies. These parameters are defined in such a 
manner that the SFR is directly related to the L TIR (Kennicutt 1998) whereas the SFE 
is the L T ir normalized by the molecular hydrogen mass (Young 1999). It has been noted 
that along the Hubble sequence, the SFR increases toward late-type, which could be due 
to dynamical instability increasing with decreasing bulge-to-disk ratio (Combes 2001). This 
trend is generally observed among the optically unbarred spirals in the present study as 
well, whereas galaxies with warmer FIR colors and hence higher SFRs are less likely to be 
transitional (SAB) or barred (SB) galaxies (See Fig.2) 

Question can also be raised such as, what types of morphological features do the cool 
and warm galaxies have and how could these features be linked with the FIR color-luminosity 
correlation ? We find that the LIRGs in our sample are mostly late type unbarred spirals. 
The transitional and barred spirals, in general, fall into the class of low luminosity SFGs. 
Among LIRGs, the late type unbarred spirals are predominately warm, and the early type 
unbarred and barred are systematically cooler. One probable reason for the transitional and 
barred spirals being low luminosity star forming galaxies is that they are in the post-starburst 
phase, where the large amount of gas driven inward by the bar have already been converted 
into stars, consistent with the findings of Sheth et al. (2005). 

Finally we note that the low luminosity early type unbarred (SA) and transitional (SAB) 
spirals are responsible for the large scatter in the C — L TIR diagram. It may be that, in 
spite of our selection process, we have dust enshrouded AGN in the sample. These sources 
could span the entire luminosity range considered and could contribute to the scatter in this 
diagram. On the other hand this scatter may also arise from more variability in the SFR 
and star formation behavior, e.g., short episodic starburst or more possible configurations, 
e.g., large bulges and small star formation disks of early type sources. 
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6. Conclusions 

We have analyzed a sample of 213 nearby IR galaxies to study the correlations of galaxy 
FIR color vs. several photometric and physical parameters. The set of measures include total 
IR luminosity, radio surface brightness as well as radio, NIR, and optical sizes. Our objective 
is to understand cool and warm IR sources using various correlations. We find that galaxy 
radio surface brightness is well correlated with FIR color whereas size is less well correlated 
with color. A weak color-size correlation signals a significant uncertainty associated with 
the distance measurements of sample galaxies. It may also reflect methodological differences 
in estimating angular size in different wavelengths. We also find that late type galaxies 
(from Sc and beyond) of all morphology classes show less scattering than early types in the 
color-luminosity diagram. Our study shows that dust in the cool IR sources is probably less 
centrally concentrated than in the more typical luminous and ultra-luminous IR galaxies 
that are dominated by more centrally concentrated warm dust. We believe this result has 
significant implication in terms of FIR color-luminosity diagram: in absence of any structural 
information of high- ,2 luminous sources one can get rough estimate of the extent of the disks 
if one simply knows the FIR colors and luminosities of galaxies. 
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Fig. 1. — FIR color as a function of L TIR for the IRAS > 1.2 Jy sample of ~ 4350 

galaxies (panel a), and our sample of 213 galaxies (panel b). Three classes of galaxies, in 
general, are shown by small triangles (SFGs), small dots (LIRGs), and stars (ULIRGs). The 
cool and warm sources are indicated, respectively, by large filled symbols and open symbols 
(triangles for SFGs and circles for LIRGs) in the regions delineated by dotted closed curves. 
In each luminosity bin the horizontal bar and (small) cross represent median and mean color. 
Each luminosity class (delineated by long vertical lines) contains equal number of cool and 
warm sources: 15+15 for SFGs and 10+10 for LIRGs. Large crosses and asterisks represent, 
respectively, local sub-mm and interacting/disturbed galaxies. A representative error bar is 
shown at the top of each panel. 
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cool systems are extended over larger area compared to the warmer one where the emission 
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Table 1. Partial list of low luminosity star forming galaxies (SFGs) 
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NGC1204™ 





.25 


1.10 


7.33 


10.18 


14.2 


58.51 




10.88 


66.0 


18.0 


33.3 


1.9 


1.0 


1.8 


SO/a: 


149 


NGC1022™ 





.71 


3.28 


19.71 


27.33 


26.4 


19.33 




10.35 


144.0 


120.0 


56.3 


4.0 


3.0 


6.0 


SB(s)a 


150 


NGC3885™ 





.57 


1.47 


11.89 


16.25 


23.4 


22.93 




10.27 


144.0 


60.0 


56.8 


8.0 


3.0 


7.0 


SAB(r:)0/a 


151 


NGC1266™ 





.25 


1.20 


13.13 


16.89 


75.7 


28.86 




10.46 


90.0 


60.0 


32.8 


5.0 


2.0 


6.0 


SB(rs)0 


152 


NGC3597"" 





.67 


2.18 


12.84 


16.21 


34.8 


48.31 




10.97 


114.0 


90.0 


25.0 


8.0 


5.0 


7.0 


S0 + : 


153 


NGC1222™ 





.50 


2.28 


13.06 


15.41 


43.8 


32.26 




10.60 


66.0 


54.0 


22.9 


11.0 


7.0 


18.0 


S0- 


154 


NGC0839™ 


(1 


.52 


2.27 


11.67 


13.03 


15.1 


51.10 




10.97 


84.0 


42.0 


29.2 


2.6 


1.3 


1.8 


SApec 



NO 



Note. — c :cool; w :warm; 0jy :major axis; m :minor axis 
Col.l :Number; Col. 2 :Galaxy Name; Col.s 3-6 :IRAS flux densities (Jy) from Revised Bright Galaxy Sample (R.BGS; Sanders ct al. 2003); 
Col. 7 :1.49 GHZ flux density (mjy) (Peak flux Sp in Condon, Anderson Sz Hclou 1991; Wang &c Helou 1992) Col. 8 : estimated galaxy 
distance d in Mpc; col. 9 : flags P.S.V on distance estimate (see RBGS Table 1 for details); Col. 10 :total IR. luminosity log-^ ot^T/i?.] 
in Lq (RBGS); Col.s 11-12:_B — band 25 mag arcscc 2 major and minor axes (arcscc); Col. 13 :K S band 20 mag arcsec - 2 major axis 
(arcsec); Col. 14-15 :1.49 GHZ major and minor axes for the resoluation at which the galaxy is resolved (arcsec); Col. 16 :1.49 GHZ beam 
rcsolutiton; Col. 17 : Galaxy optical morphology. 
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Table 2. Correlation statistics from Pearson's Correlation Test 



Galaxy Class 


C- 


-sb r 


C 


-D R 


C - 


- Dnir 


C 


-D B 




r 


p 


r 


P 


r 


P 


r 


P 


SFGs+LIRGs+ULIRGs 
SFGs+LIRGs 


0.77 
0.77 


1.0E-08 
1.0E-08 


-0.65 
-0.67 


2.14E-07 
1.83E-08 


-0.14 
-0.22 


4.74E-02 
1.91E-03 


-0.11 
-0.19 


1.15E-01 
5.33E-03 



Note. — Correlation statistics for FIR color C vs. all other parameters in Fig. 2. In 
this table, r is the correlation coefficient, and P is the probability that correlation could 
arise from an uncorrelated sample. The confidance is 1 — P. Galaxy 1.49 GHz surface 
brightness and size are represented by SBr and Dr. For both K s and B bands, test is 
performed after removing galaxies with unknown sizes. 



